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We review the difference between standard environmental decoherence and 'intrinsic decoherence', 
which is taken to be an ineluctable process of Nature. Environmental decoherence is typically 
modeled by spin bath or oscillator modes - we review some of the unanswered questions not captured 
by these models, and also the application of them to experiments. Finally, a sketch is given of a 
new theoretical approach to intrinsic decoherence, and this scheme is applied to the discussion of 
gravitational decoherence. 
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I. INTRODUCTION 



In this short paper I was given the job of briefly review- 
ing current progress on the theory of environmental deeo- 
herence as applied to condensed matter systems. In order 
to make the topic a little more interesting I have added 
some material on 'intrinsic decoherence'. This latter term 
denotes a process which is considered to be generic to Na- 
ture, and which causes phase decoherence even without 
any 'averaging over an environment'. As such it is equiv- 
alent to a breakdown of quantum mechanics. Many argu- 
ments have been given to suggest that such a breakdown 
is something to be looked for - perhaps the most com- 
pelling are those beginning either from the clash between 
quantum mechanics and General Relativity, or from the 
difficulties we have in thinking about macroscopic super- 
positions. 

We will see that the topic of intrinsic decoherence is 
in its infancy, and that it is then interesting to compare 
it to environmental decoherence where things are a lot 
clearer. However in what follows I will also try to em- 
phasize some of the remaining problems (as well as some 
of the successes) , that we have with environmental deco- 
herence. 



II. ENVIRONMENTAL DECOHERENCE: 
EXISTING MODELS 

In order to say something useful about environmental 
decoherence, let's begin by reviewing the current under- 
standing of how this works. This has come a very long 
way since the schematic models that were used before the 
1980's. One begins by writing the universe as a Hamil- 
tonian of form 
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Ho{P, Q)+Hint{P, Q; Pi,,qu)+Henv{Py, Qu) (1) 



where the Q, P are generalized coordinates and momenta 
for the central system of interest, and the {qu,Pu} are the 
same for the environmental degrees of freedom. The vari- 
ety of models for environmental decoherence then comes 
from: 



(i) Models of the Environment: At first glance it 
seems as though we ought to require a large variety of 
different model Hamiltonians for different environments. 
However with few exceptions, it has been possible to re- 
duce all of these to one of two models, as follows: 

(a) Oscillator Baths: Here we model the environment 
by a bath of independent oscillators [H 0] , so that 
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(2) 



where the {tOq} are the frequencies of the No oscillators. 
Opinion on oscillator bath models was polarized in the 
early days - sometimes they were applied willy-nilly to 
almost any environment. And yet in their original dis- 
cussion of such models, Feynman and Vernon [l| made 
it quite clear that the oscillator bath representation was 
only generally valid if the coupling between each environ- 
mental mode and the central system was very weak. The 
most obvious case is when the environmental modes are 
delocalized, extending over some large volume (so that 
No is an extensive quantity, depending in general on the 
UV cutoff r^o we employ). Examples are phonons, pho- 
tons, magnons, electron-hole pairs, and other collective 
excitations in a fermionic system. The coupling is then 
taken to have the general form P, Q 



TTOSC 



N 
9=1 



[Uq{Q)Xq +Vq{P)pq 



(3) 



which is restricted to being linear in the oscillator coordi- 
nates Xq and their conjugate momenta Pq. The weakness 
of the interactions Uq{Q), Vq{P) is assured by the delocal- 
ization of the modes (so that each mode has a normal- 
ization factor ^ No ^ , and thus \xq\, \pq\ ~ m 
the 'thermodynamic limit' (ie., for large systems), when 
A^o ^ 1, this works. 

Actually this model is far more generally applicable 
than it seems. Often couplings to oscillator modes that 
are higher-order in the oscillator coordinates are impor- 
tant or even dominant in the real world - but these can 
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be lumped into linear couplings to effective oscillators 
via temperature-dependent couplings. And while anhar- 
monic couplings between oscillators will obviously affect 
the very long-time dynamics of the bath, they typically 
have little effect on the system dynamics over experimen- 
tally relevant timescales - which is what we are interested 
in. The basic idea behind the oscillator bath model is one 
of linear response by the bath modes, to the perturba- 
tion by the system - as such the idea goes back at least 
to Rayleigh's model of friction 01 . 

(b) Spin Baths: Where oscillator bath models typically 
fail is when the environmental modes are not delocal- 
izcd. This is quite common - examples include lattice de- 
fects, dislocation modes, spin impurities, dangling bonds, 
nuclear spins, and localized phonons/vibrons; there arc 
many others. Because these modes are localized in space, 
their Hilbert space is very small - most commonly they 
can be modeled as a spin bath' of 2-level systems [5[ , de- 
scribed by a set of Ng 'pseudospins' {cTk}, with a Hamil- 
tonian 

AT, AT, 

iJlL = Y.^k-<Tk+ Y.'^^f.atal (4) 

where each pseudospin feels a set of external fields {hfc}, 
and inter-spin interactions Vuk' ■ These arc then coupled 
to the central system by an interaction: 

HZ = ^F,iP,Q)-a, (5) 

k 

Unlike the case of the oscillator bath, the interactions 
Fk{P, Q) are not necessarily weak; there is no reason to 
assume that \Fk{P,Q)\ N7 ' (indeed m many cases, 
Fk{P, Q) is independent of -/V^Oj ^"^^ moreover, the local- 
ized modes may also be limited to some region in space, 
so that their number Ng is limited. We note also that 
one often must include the interactions V^^, between the 
bath spins; even though they are weak, they can have 
significant effects if Ng is not large. A key dimensionless 
parameter for a spin bath is Afc = |i^/c(P, (5)|/|hfe|, which 
tells us the relative strengths of the two fields acting on 
the bath spin a^. If Afe 3> 1, the bath spin dynamics 
is slaved to the central system; in many experimental 
systems this is the case. 

We see that the spin bath is inherently quantum- 
mechanical: discrete environmental states in a finite 
Hilbert space have no classical analogue. At low tem- 
peratures, such modes are everywhere in a solid - even 
very pure and well-prepared metallic or semiconducting 
crystal have a non-negligible concentration of defects and 
dislocations, as well as paramagnetic impurities, and typ- 
ically their surfaces are very hard to control and will have 
surface impurities, dangling bonds, and various kinds of 
lattice defect and 'charge fluctuator' modes, all of which 



can be described as spin bath modes. And this of course 
is quite apart from nuclear spins, which exist in profusion 
in almost all solids, and local vibrational modes, which in 
many non-crystalline and organic systems play a key role, 
coupling strongly to, eg., electronic modes, and behaving 
at low T very much like 2-level systems. 

Reviews of both oscillator bath and spin bath models 
can be found in the literature 0, Q . 

(ii) Models for central systems: These of course 
abound, and theory has examined a large number of 
them, both for their intrinsic interest, and in the con- 
text of experiments. The best known class of models is 
one in which the central system is reduced to a 'qubit', 
ie., a 2-level system with Hamiltonian 

Ho = Bo(t) • T (6) 

where Bo(<) is some 'external' time-dependent field act- 
ing on the qubit. Most theoretical analyses of decoher- 
ence for this system have been done assuming a static 
Bq. When coupled to an oscillator bath, via a cou- 
pling = '^^UqxA2,_we have the diagonally-coupled 
'Spin-boson' model [a, LZI; it has been studied in the 
weak-coupling limit ever since the beginning of quan- 
tum mechanics (notably in atomic physics and quantum 
optics), and for more general couplings starting in the 
late 1960's, in the context of RG theory and the Kondo 
problem. One can also add a non-diagonal coupling 
H^it = Tliq^q^+H-c^Xq^ which cau be important. 

When coupled to a spin bath, via a coupling = 
F/c • (TkTzi the coupled system is known as the 'Cen- 
tral Spin' model Such models were first studied in 
NMR and ESR, where again work was limited to weak 
coupling; more general studies began with work on tun- 
neling magnetic systems [1], and have accelerated with 
work on decoherence in solid-state qubit systems [^, 0] . 
Again, one can generalize this to include non-diagonal 
terms, leading to a general tensor interaction. 

The Spin-boson and Central Spin models have very 
different dynamics 0] ; this is inevitable, given the en- 
tirely different structure of the Hilbert spaces for the 
two baths. The differences are seen most easily in the 
lineshape (eg., in x"z(w), the Fourier transform to fre- 
quency space of the time correlation {Tz{t)fz{Q)))- Sup- 
pose we hold Bo(i) constant in time. Without the baths, 
x"z(a;) oc 5{ijj — Wo), where ujo = ^Bo- As one increases 
the coupling to an oscillator bath, the ^-function for the 
spin-boson model first spreads into a narrow Lorentzian, 
and then with increased coupling, the peak broadens and 
moves to lower frequency; if the coupling becomes strong 
enough (as in an Ohmic bath at strong coupling) the peak 
moves to zero energy (critical damping), and eventually 
collapses to a (5- function at zero frequency (ie., localiza- 
tion in time). On the other hand for a spin bath, the 
Central Spin lineshape typically becomes very unconven- 
tional as one increases the coupling, developing multiple 



3 




0.0 0.5 1.0 cO„ / A„ 

FIG. 1: Plot of 7m Xzz{i^), the absorption probability de- 
scribed in the text, for a qubit coupled to a spin bath (the 
'Central Spin' model). Without decoherence this would be a 
sharp line at lj/ Ao = 1. The shaded part shows contributions 
from processes where the spin bath polarisation is unchanged. 



peaks and looking anything but Lorentzian. The cor- 
responding behaviour in time also looks very different. 
Fig. [T] illustrates this for the Central spin model in a 
case where the dimensionless coupling to the spin bath is 
intermediate in strength - the contributions are divided 
between processes in which the net bath polarization does 
not change, and those for which it does. For more de- 
tailed discussion of such results see ref. Q . 

The reason for this difference is that (i) spins have a 
dynamics very different from that of oscillators (the spin 
dynamics is highly constrained compared to that of the 
oscillators), and (ii) in the Spin-boson model the bath is 
always in the linear response regime - each oscillator is 
only weakly affected by the qubit - whereas in the Central 
spin model, the bath spin dynamics is strongly affected 
by the qubit, and we are typically far from from a linear 
response regime. If Afe 3> 1, the bath spins simply follow 
the central qubit (with deviations ~ 0(A^^)). 

Another key difference between the Spin-Boson and 
the Central spin model is related to this latter feature. 
In the spin-boson model (and indeed in all oscillator bath 
models) dissipation and decoherence are closely linked 
0; one cannot have one without the other, and both 
are connected to bath fluctuations via the fluctuation- 
dissipation theorem. However in the Central spin model 
it is not only possible but quite typical for strong deco- 
herence to occur with little or no dissipation. Indeed, the 
main decoherence mechanism in the Central spin model 
is that of 'processional decoherence', in which the bath 
spins process in the time- varying field of the qubit, with- 
out any energy exchange between them 0, [l^ . Again, 
this is a fairly generic feature of spin bath decoherence. 

For a more extended discussion of the differences be- 



tween spin and oscillator baths, see ref. [lO[. 

Many other models for central systems have been dis- 
cussed in the literature, including tunneling systems, 
multiple qubit systems, quantum critical systems, and 
quantum information processing systems. We will only 
look at one, viz., a set of coupled qubits, with Hamilto- 
nian 



(7) 



in which each local field Bj(i) on the j-th qubit is sub- 
ject to independent control, as are the couplings J'^f {t) 
between them. Then, to model decoherence, we couple 
each individual qubit to the oscillator bath and the spin 
bath, using the same local couplings as given above for 
a single qubit (but with, in general, different couplings 
on each qubit). There are many possible variants here - 
different qubits may couple to different spin baths (being 
localized, these baths may be quasi-independent), and we 
by changing the form of the coupling to oscillator baths 
(Ohmic, super-Ohmic, etc.) we can radically alter the 
dynamics. 

However a new kind of 'internal decoherence' enters 
with the introduction of the inter-qubit interactions. 
Much depends on the form of the Jij (in particular, 
whether they are long-range or short-range in space) ; we 
will be interested particularly in long-range dipolar inter- 



actions, which we can write in the form J 



a/3 



j: 



0Tn,a/3 



where the interaction strength is J?- = ^ i^'s I foi' 
spins separated by a distance r^j, and the spin tensor 
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where the are the effective "g- factor" tensors for the 



qubits [25|. Now these interactions are 'frustrating' [11 



ie., they make it hard for the system to order, and hard 
for simple pairwise interactions to operate locally; and 
they are long-range. The radical effect that these inter- 
actions can have is illustrated in Fig. [31 which shows 
a set of typical 'resonating pairs' of qubits, able to per- 
form 'flip-flop' transitions via the interaction J^. The 
key point here is that in a 3-dimensional system, such 
pairs are just as likely to be found far apart from each 
other as nearby (the decrease in interactions strength be- 
ing compensated by the larger number of available pairs 
at longer range) . Thus the quantum dynamics of the sys- 
tem is not driven by short-range 'quantum diffusion' of 
information - the information diffusion is spread all over 
the system. 

Why do we care about this? One reason is that it 
is very hard to get rid of dipolar interactions like this - 
they are not easily screened. So quantum computation is 
faced here with the problem that unless one can actually 
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FIG. 2: (Color on-line) The dipolar interactions between 
qubits can cause resonant flip-flop transitions over long ranges 
- the figure depicts the qubits for which such a resonant in- 
teraction might exist in some qubit network 



use these interactions as part of the computation pro- 
tocol itself, then they cause 'correlated errors' over long 
ranges. This is of course just another form of decoher- 
ence. Thus we see that Hamiltonians like ([7]) containing 
dipolar interactions describe a key aspect of all qubit- 
based quantum information processing systems. 

Study of problems like this on real systems has until 
recently been limited to looking at the effect of both the 
dipolar interactions and the coupling to the spin bath 
on, eg., the quantum phase transition in these systems 



(see, eg., ref. 



12|; note that relat ed p roblems have been 
discussed in the NMR literature [l3|). However in sec- 
tion IV we will discuss how it has now become possible 
to quantiiy the various decoherence processes in such a 
system. 



III. ENVIRONMENTAL DECOHERENCE: SOME 
QUESTIONS 

A number of questions about the preceding framework 
jump out at us. here are some: 

(a) How good are these models? Thus - under what 
circumstances are they valid, and how does one derive 
them for given systems? And when can one use an oscil- 
lator bath representation even for localized modes? What 
about models like 'noise' models, where we model the en- 
vironment by a noise source? And so on. 

(b) Are there other kinds of decoherence (not captured 
by these models)? What about decoherence caused by 
time variation of parameters, or decoherence existing in 
the interaction channels - neither of these is captured in 
the above. What about 'intrinsic decoherence, ie., deco- 
herence quite independent of any environment, intrinsic 
to Nature? 



(c) What about Experiment? Can we have large- 
scale quantum computing, or large-scale quantum de- 
vices? What are the real limits on quantum behaviour 
imposed by decoherence? And - is decoherence inevitable 
at the macroscale? 

(d) What about the 'Quantum Measurement Prob- 
lem' ? If there is such a problem, does decoherence 'solve' 
it, or any of the other foundational problems in quantum 
mechanics? 

In the rest of this section we take a whirlwind ride 
through questions (a) and (b) . The next section will look 
at (c); and the rest of the paper will take a closer look 
at intrinsic decoherence. As for question (d), my general 
feeling is that until questions (a)-(c) have been resolved, 
there is really nothing new that can be said about (d). 



(a) How good are these Models?: The derivation of the 
oscillator bath models was discussed by Leggett et al. 0, 
|3| . The key here is the assumption of adiabaticity - if we 
define instantaneous eigenstates (f)n{x,y,Q) and energies 
^n{Q) for the environment when the central system has 
coordinate Q, then the oscillator energies coq = cjnm = 
^niQ) — ^miQ), and wc can define the couplings Uq, Vq in 
terms of the gauge potential 

n 



dx,. 



d ~ 

^{Xu)-g^(|>rn{Xu) (9) 



provided that Q{t) varies slowly enough in time that 
l^nml ^ \{^n — ^m)\- If this couditiou is obeyed, then 
we can use the oscillator model even if the original cou- 
plings to the bath are strong - the {uq,Vq} will still be 
small. 

However in many cases there will be low-energy envi- 
ronmental modes for which this condition is not obeyed 
(indeed some of them will 'resonate' with the central sys- 
tem, ie., their frequencies will match the characteristic 
frequencies of dynamics of the central system). In this 
case one must fall back on the perturbative assumption 
discussed above, viz., that the the actual couplings to the 
bath modes are weak. It is precisely this assumption that 
fails in the case of many localized modes. 

How good then are the spin bath models? In some 
cases the couplings Fk{P, Q) and fields hfc are known ac- 
curately from microscopic arguments and/or experiment 
(eg., with hyperfine couplings to nuclear spins, in spin 
qubit systems). In other cases we can sometimes only 
guess at these quantities (as with the coupling to defects 
in magnets or superconductors). The key assumption 
however is that the couplings Fk{P, Q), or their rates of 
change, are small enough so that the bath modes are not 
excited outside their restricted Hilbert spaces. Note, inci- 
dentally, that we also assume that V^f!! is small compared 
to either |hfc| or \Fk{P, Q)\; otherwise the bath spins cou- 
ple together into delocalized modes, which behave more 
like oscillators. 

It has quite common in some fields to approximate the 
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effect of an environment by coupling the system to an ex- 
ternal 'noise' source, ie., to approximate the full Hamil- 
tonian in (P by a form H{t) ^ Ho{P,Q) + VN{Q;t), 
where the time-dependent potential V^r contains a ran- 
dom force acting on Q, and the effects of this force are av- 
eraged over, using some technique (Langevin, Lindblad, 
Fokker- Planck, etc.). But it well known that such mod- 
els can never fully capture the effect of integrating out a 
dynamical environment . They only describe the real 
part of the influence functional, which itself completely 
characterizes the effect of the dynamic environment on 
the system [l|,[l3; moreover, when approximated by, eg., 
Markovian processes, they also lose the long-time correla- 
tions which are characteristic of the influence functional 
of many environments. Physically, it is clear that a key 
part of the environmental effect on the central system 
dynamics is via a time-retarded 'back-reaction' on this 
system. We note that the results may sometimes be quite 
counter-intuitive. The example of a particle moving on 
a lattice while coupled to a spin bath provides an 
example: the particle density matrix shows simultaneous 
ballistic motion and anomalous diffusion, reminiscent of 
weak localization. Further discussion of all this is quite 
technical. 

A related question is that of the crossover between spin 
bath and oscillator bath models as one switches on the 
coupling between the bath spins - when the interspin in- 
teraction dominates, we expect the spin bath to behave 
like an oscillator bath. This crossover (which for a macro- 
scopic bath should under some circumstances become a 
phase transition) has been discussed in some detail for 
the Central Spin model but there has been no gen- 
eral analysis as yet. 

(b) What other kinds of Decoherence exist?: A more 
obvious question about failures of the oscillator and spin 
bath models is - what do they leave out? This turns out 
to be an interesting line of investigation. Here are a few 
suggestions (the list is not meant to be exhaustive): 

(i) Other environmental modes: It is certainly not 
clear that the set of all environmental modes that might 
cause decoherence is exhausted by oscillators and local- 
ized 'spin bath' modes. One possible example, found in 
many systems at low energies, is that of of delocalized 
quantum soliton modes - it is by no means clear (and 
has certainly never been shown) that these can be de- 
scribed in all cases by an oscillator bath model. There 
may also be other possibilities 

(ii) 3rd-party Decoherence: Since environmental deco- 
herence arises from the entanglement of a system with its 
environment, we can certainly envisage situations where 
this happens without the direct interaction of the two. 
For example, the state of each may be 'prepared' or oth- 
erwise conditioned by a 3rd party. A simple example 
is that of the coupling between the vibrational and ro- 
tational modes of a buckyball, on the one hand, and its 



centre of mass motion, on the other, via the slits it passes 
through in a 2-slit experiment [lo| . But one can easily 
imagine more subtle examples, and the real question of 
interest here seems to be - how can we know that a sys- 
tem has not previously entangled with some environment, 
perhaps long in the past, perhaps via a 3rd party, in such 
a way as to cause decoherence now? This seems to be to 
be an open question, relevant in principle to many ex- 
periments. 

(iii) Switching Decoherence: It is quite common, par- 
ticularly in quantum information theory, to write time- 
dependent Hamiltonians - we have done so above in (|6]) 
and ([7]). But this raises two points. First, where is the 
time-dependent field coming from that is causing this 
time variation, and shouldn't this source term also be 
included in the Hamiltonian (to take account of possible 
back- reaction in the source)? And second - how do we 
know that the time variation cannot project the system 
dynamics out of the Hilbert space of the effective Hamil- 
tonian we are using? This is of course the whole problem 
in 'adiabatic quantum computation', and in the similar 
problem of the dynamics of quantum phase transitions, 
where a large number of energy gaps can decrease to very 
small values during the time evolution. In this case, one 
can also include all these extra levels in the 'environment' 
(now with time-dependent energies). This topic is ripe 
for investigation, and certainly there are many papers on 
the dynamics of driven systems coupled to a bath, in- 
cluding some discussions of decoherence for simple mod- 
els [l6[. However, to my knowledge no systematic study 
is yet available. 

(iv) Channel Decoherence: in quantum information 
theory it has been common to discuss the decoherence 
arising during the transmission of information down a 
channel. But such decoherence is ignored in most of the 
effective Hamiltonians we write down, where coupling 
constants are just treated as constants - this is equiv- 
alent to assuming instantaneous interactions. The value 
of these interaction constants depends on the UV cutoff 
hrio we assume in the theory - but this takes no account, 
for a finite quantum computer, of the time To taken for 
signals to cross the computer. In reality there is a set 
of timescales - the time to taken for qubits to switch, 
the time To = I/^Iq which is a UV cutoff time, and the 
timescales associated with the qubit operations them- 
selves. Unless To is very short compared to all of these, 
we should not be assuming an instantaneous interaction, 
but rather a retarded propagator. A proper treatment 
of this will be equivalent to a microscopic derivation of 
channel decoherence. 

(v) Intrinsic Decoherence: By intrinsic decoherence I 
will refer to a decoherence in Nature that is ineluctable; 
it does not result from averaging over any kind of ex- 
ternal environment, and thus can even exist for an iso- 
lated system. From this we see that it must amount 
to a breakdown of quantum mechanics, because the lin- 



6 



ear equation of motion on quantum mechanics preserves 
superpositions for all time. The possible breakdown of 
quantum mechanics has been investigated in many dif- 
ferent ways in the last 40 years; some of the best known 
efforts include the postulation of additional non-linear 
terms 17, 18 1, or else stochastic terms 0, [1^, to the 
Schrodinger equation, and a breakdown of this equation 
caused by gravitational effects [2l|, 122] ; this list is by no 
means exhaustive. In most cases the effect of these devi- 
ations from standard quantum mechanics is to mimic de- 
coherence - in other words, in an experiment they would 
be expected to add a further contribution to the existing 
environmental decoherence, which may be very hard to 
distinguish from it. 

It is fair to say that our understanding of all five of the 
above decoherence mechanisms is more or less open right 
now - and all of them are interesting and potentially im- 
portant for experiments. In the second half of this paper 
I will be looking at just one of them, viz., intrinsic de- 
coherence. However, a key challenge for any experiment 
looking for intrinsic decoherence will be to distinguish 
all sources of environmental decoherence from any intrin- 
sic decoherence (and if possible eliminate these sources) . 
Therefore we first take a look at how it is possible to nail 
down environmental decoherence sources in experiments 
(cf. question (c) above). 



IV. ENVIRONMENTAL DECOHERENCE: 
THEORY VS. EXPERIMENT 

In the last decade many experiments have been done 
looking at decoherence in both quantum optics and con- 
densed matter systems - in the last couple of years these 
have been extended to hybrid systems like optomechan- 
ical, optomagnetic, or opto-superconducting systems, in 
which some combination of optical cavities, supercon- 
ducting or magnetic qubits, and mechanical or electrome- 
chanical oscillators are coupled and put collectively into 
a quantum state. In some of these experiments the de- 
coherence mechanisms were fairly simple, involving cou- 
pling of an oscillator degree of freedom to photons or 
phonons, and in this case the experimental decoherence 
rates did not wildly disagree with simple theory. In more 
complicated cases, usually involving larger systems where 
the environment involved spin bath degrees of freedom, 
the experimental decoherence rates have often been much 
larger than naive theory would predict. 

The temptation is to simply write these discrepancies 
off as 'junk' effects, coming from some combination of de- 
fects, impurities, etc. However the problem here is that if 
the experiments are supposed to be testing, eg., the appli- 
cability of quantum mechanics to large and/or complex 
systems, such a temptation is no longer open to us - the 
theory should then also be predicting the environmen- 
tal decoherence rates, and a genuine experimental test of 



quantum mechanics will then be testing such predictions. 

It is of course always possible to argue that the the- 
oretical models are perfectly OK, but that we are an- 
alyzing their dynamics incorrectly. All theoretical pre- 
dictions assume that we can average over all environ- 
mental modes to get a reduced density matrix for cen- 
tral system alone; the dynamics of this reduced den- 
sity matrix is then described by a propagator K(2, 1) = 
K{Q2, Q'2] Qi,Qi]t2, ti), so that the reduced density ma- 
trix p{2) = p{Q2, Q'2] h) evolves according to 



P(2) 



dlif(2,l) p(l) 



(10) 



where dl means / dQidQ'^. The propagator can be writ- 
ten in the general path integral form [l| 

A'(2,l) = VQ /^'PQ' e^(^°l^l^^°[^'l)j-[0,g'] 

(11) 

where So[Q] is the free central system action, and 
J^[Q,Q'] is the "influence functional", defined in general 
by 



nQ,Q']^l[(^AQ,t)uliQ',t)} , 



(12) 



Here the unitary operator U^{Q,t) describes the evo- 
lution of the i/-th environmental mode, given that the 
central system follows the path Q{t) on its "outward" 
voyage, and Q'{t) on its "return" voyage; and T[Q,Q'] 
acts as a weighting function, over different possible paths 
iQit),Q'{t')). ^ ^ 

The calculation of the influence functional is not nec- 
essarily easy. It is given in general form by a functional 
integral over the environmental modes, as 

=n / T^qlc^p'-^iq.,ql;Q,Q') (13) 

where the complex phase ^'(9,^, ql; Q, Q') is given by the 
difference in actions over the two paths in the reduced 
density matrix, ie.. 



(14) 



The 'bath action' S's ((7^,(5) just sums over the environ- 
mental and interaction terms, ie., we have Sb = / drLs, 
where 

Lsiqi^, Q) = Lint{qu, Q) + Lenv{{q^) (15) 

Everything then depends on the form taken by the 
influence functional. The oscillator bath influence func- 
tional was derived by Fcynman and co-authors [l|, Q ; one 
lets q^ Xq, and the Lagrangian Lsixq, Q) is then 



Lq {Xq , Xq , t) 



''^±1 -UJqX^] -Jqit)Xq 



(16) 
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where jqit) is the driving force acting on the q-th bath 
oscillator, coming from the central system. Thus, eg., 
if the coupling in ^ is Li^t = ^J2q^giQ)-''qj ^® have 
7g(t) = Fq[Q{t)). The bath action 5*^ is then a quadratic 
form in the oscillator coordinates and one can find an 
exact expression for iF[Q,Q']. However, depending on 
the form of Fq(q) as a function of Q, the final integrals in 
(fTTj) may still be highly non-trivial. This topic has been 
studies at great length - a good review is by Weiss . 

The spin bath influence functional is very different, and 
has been studied much less. Now one lets ^ a^, and 
the Lagrangian LsicTk, Q) takes the form 



L{cTk,CTk\t) = Ak 



dT 



Tk{t)-ak (17) 



where Ak is formally the static gauge field which is pro- 
duced by a monopole of unit strength at the centre of the 
Bloch sphere for cr^ (and which leads to the Berry phase 
for (Tk |23|). The time-dependent field Tk{t) is given by 



(18) 



in which hfe is the static external field acting on the k-th 
bath spin, and Fk{t) is the dynamic field acting on this 
spin, coming from from the central system (c/. eqtns. ([4]) 
and ([5])). The extra i-dependent fluctuating field ^^(t) is 
defined by its components as 



it) 



(19) 



and to first approximation this can be dropped, since 
V^i^ is a small perturbation on the field hj. -I- Ffc(t). 

This Lagrangian is much harder to deal with than the 
oscillator Lagrangian - indeed the solution to the prob- 
lem of a spin- 1/2 in almost any but the simplest time- 
dependent fields is practically impossible to determine 
(something which always surprises students). In spite 
of this, detailed expressions have been produced for the 
dynamics of systems coupled to spin baths, which are be- 
lieved to be accurate in a fairly wide parameter range, at 
least up to a certain time 0, Nevertheless, given the 
uncertainties here, one obvious way that one might try 
to explain the discrepancies between theoretical and ex- 
perimental environmental decoherence rates in so many 
experiments, is simply to argue that the theoretical cal- 
culations of spin bath decoherence rates are just way off 
target. 

What is obviously required here is a confrontation be- 
tween theory and experiment where one has accurate 
knowledge of the physical sources of the decoherence (no- 
tably the spin bath modes). In what follows I would like 
to briefiy describe one specific experiment, which did ex- 
plicitly attempt such a test. The experiment looked at 
a macroscopic crystal of tunneling magnetic molecules, 
and specifically was looking for the macroscopic spin po- 
larization wave that should propagate in a region of the 



sample if the decoherence rate is low. The experiment 
is interesting because it checks out not just decoherence 
from local coupling to both oscillator and spin baths, but 
also a more insidious decoherence which comes from the 
long-range dipolar coupling between all of the tunneling 
molecular spins. 

The experiment was done on a crystal of "Fe-8" 
molecules; each molecule contains 8 Fe ions in a small 
core, which are surrounded by a complicated sheath of 
ligand molecules, containing H, C, N, O, and Br ions 
(see Fig. [3]). The spins in the central magnetic core 
are locked together by superexchange interactions into 
a spin-10 system, whose dynamics is governed by a bi- 
axial crystal field. Depending on the different isotopic 
concentrations in the molecules, each one contains up 
to 216 nuclear spins. For the discussion here we can 
just imagine that at sufficiently low temperature (below 
about 2.5 K), we have a set of spin qubits, each of which 
couples locally to both a phonon oscillator bath and a 
nuclear spin bath. These couplings, including all of the 
hyperfine couplings to the different nuclear spins, are very 
well understood quantitatively. The same is true of the 
dipolar coupling between the qubits, and so we can model 
the system very accurately by the effective Hamiltonian 
H = Ho + HsB + Hose, where 



(H^)^° 



(20) 



sums over all crystal lattice sites, and where Hsb and 
Hose describe the nuclear spin and phonon terms respec- 
tively. The temperature T of the system is defined by the 
phonons; there is also an magnetic field Hj^ applied in 
such a direction that we can control the 'tunneling split- 
ting' Ao over a very wide range (some 7 orders of magni- 
tude), without adding any terms cx fj. The effective in- 
teraction energy J"^^ (H.±)t"tj is actually dominated by 
longitudinal terms Jfj^, so Ho then becomes the famous 
dipolar quantum Ising model. If the interaction were 
short-ranged, with strength Jo, we would get a standard 
zero-T quantum phase transition when \Ao/Jo\ ^ 0{1). 
However the interaction here is dipolar, and moreover we 
also have the couplings to the oscillator and spin baths - 
this changes things fundamentally. 

This system is a key test case for quantum computing, 
for the following reason. We expect almost any solid- 
state gate-based quantum computer to possess local cou- 
plings between the qubits and spin and oscillator baths; 
and moreover, we expect long-range dipolar interactions 
to exist between the qubits, either electric or magnetic (or 
even strain-mediated). Such long-range interactions are 
very dangerous for quantum computation, because con- 
ventional error correction routines do not work for them 
- indeed it is not yet clear how one would correct errors 
arising from them, although there are arguments that 
concatenated error correction could do the job j24| . In 
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FIG. 3: Simplified ball-and-stick representation of the 
[Fe802(OH)12(tacn)6]8+ system (the "Fe-8" molecule), 
showing the octahedral iron ions, the oxo and hydroxo bridges 
and the 1,4,7-triazacyclononane (tacn) ligands completing the 
coordination spheres. 
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any case, we see that this model contains all of the agents 
of decoherence that wc may expect to arise in a multi- 
qubit quantum computer or quantum register (apart per- 
haps from the switching decoherence and channel deco- 
herence noted above) . In this sense the experiment to be 
described is the first one to test all of these decoherence 
mechanisms together. A key feature of this experiment 
was that all the hyperfine interactions to the nuclear spin 
bath were known accurately, and so one can in principle 
directly test the theory of spin bath decoherence, and the 
calculational methods used to derive decoherence rates. 

The theory for the decoherence rate as a function of 
temperature, applied field, and for different isotopic con- 
centrations of different nuclear spins, was actually car- 
ried out in 2006, and gave detailed predictions for the 
Fe-8 system [iH, in a spin echo configuration suggested 
in that paper. The basic idea was that one could mea- 
sure the '1-spin' decoherence rate by doing a global spin 
echo experiment, in which the entire spin system was set 
into a spatially uniform precessional motion. One would 
of course like to measure the dynamics of multi-spin den- 
sity matrices, and this will be a key goal for the future. 
In any case, one finds theoretically that the dimension- 
less decoherence rate (defined as = l/r2Ao, where Ti 
is the 1-spin decoherence time extracted from spin echo 
measurements) should vary enormously depending on the 
field (which controls Aq), and on the temperature (see 
Fig. 21). The phonon decoherence rates increase quickly 
as the qubit tunneling frequency Aq increases (because 
the phase space for phonon emission is increasing). On 
the other hand the nuclear spin bath decoherence rates 
were predicted to fall rapidly at high Aq because the nu- 
clear spins can no longer follow the qubit dynamics. Thus 
for a single qubit, one arrives at an optimal Aq (the 'co- 
herence window' (26|). where the decoherence will be at 



FIG. 4: (Color on-line) Dimensionless spin decoherence rates 
as a function of the qubit splitting 2Ao for a crystal of Fe-8 
molecules, with a field H_L||i: along the x-axis. Brown dashed 
line - nuclear contribution to decoherence rate (upper line - 
protons at H sites; lower line - deuterons at H sites). Solid 
lines - magnon contributions at different temperatures, caused 
by intermolecular dipolar interactions. Long dashed lines - 
phonon contributions corresponding to the same tempera- 
tures as in the magnon case. Inset: the splitting 2Ao as a 
function of applied transverse field 

a minimum. However the dipolar decoherence, mediated 
by magnons, messes this up completely - and the only 
way to suppress it is by going to low temperatures (the 
qubit flip-flop processes which excite the magnons, and 
which are caused by dipolar interactions, are gapped in 
energy; they are then exponentially suppressed at low T). 

A number of attempts were made by different experi- 
mental groups to test this theory, but they were rendered 
very difficult by the strong dipolar decoherence, which 
forced a very low-T experiment when conventional fields 
were used. This problem was circumvented by Taka- 
hashi et al. [i^l, who went to fields in excess of 11 T, 
and also changed the field configuration. A very satis- 
factory agreement with theory was then found, with no 
adjustable parameters in either the theory or the experi- 
ment. This experiment actually accomplished some other 
notable goals - it was the first experiment to see macro- 
scopic spin precession in an ordered array of qubits, and 
also achieved a lower decoherence rate than any previ- 
ous molecular spin qubit system (the decoherence "Q- 
factor" or inverse decoherence rate, went to 1.5 x 10^, 
and if these experiments had gone to lower T, it should 
have gone to 6 x 10*"). However what interests us here is 
that it tested all three decoherence mechanisms simulta- 
neously (including the first ever measurement of dipolar 
decoherence). This is because by varying T and Hj^, 
one can independently vary the three mechanisms). So 
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the obvious question is - what was being tested in this 
experiment? 

The two most obvious conclusion from the results are 
(i) that since apparently all three decoherence mecha- 
nisms checked out, there could have been only negligible 
'extrinsic' effects in the experiment - no hidden decoher- 
ence from dirt effects like hidden defects, dangling bonds, 
etc., so that the only spin bath modes were the nuclear 
spins; and (ii) that the calculations of both spin bath de- 
coherence and dipolar decoherence are actually reliable. 
To this we can add that the agreement with experiment 
for phonon decoherence was to be expected - this is a 
well-understood process. The success of the predictions 
for spin bath and dipolar decoherence is welcome - for 
a long time it was not obvious if disagreements between 
theory and experiment came from experimental dirt or 
inadequate theory. 

But let us now return to the question (c) posed in sec- 
tion 3. What can we now say about decoherence for large- 
scale quantum processes, including the kind that would 
be involved in quantum computing with many qubits? 
Indeed, to what extent can one say that in doing exper- 
iments like this we are testing fundamental theories of 
quantum mechanics on large systems, in the same way 
that experiments on, say, coherence in superconducting 
SQUIDs [il] are supposed to be doing? 

It seems to me that there are two basic approaches 
that one can take to this question. One is to probe the 
predictions of quantum mechanics in situations where its 
results arc very counter-intuitive, and apparently violate 
deep prejudices we have about Nature. Examples of this 
approach are the inequalities derived by Bell [29[ , which 
probe non-local entanglement, and have now been tested 
in many different ways on microscopic systems (and in- 
deed constitute a 'resource' for quantum computation); 
and the generalization of these to time inequalities, to 
discuss possible experiments testing 'macroscopic real- 
ism' on superconducting SQUIDs [30[. It is clear that 
the experiments described above do not yet test quan- 
tum mechanics at the macroscopic level in this sort of 
way, but they certainly could be adapted to do so, and 
probably will be in the near future. 

However a second way of testing quantum mechan- 
ics, rather than just probing its consequences, is to de- 
velop alternative theories, and then do experimental tests 
which compare these directly with quantum mechanics, 
and indeed distinguish explicitly between them. This is 
precisely what theories of intrinsic decoherence aim to 
do, and so we now turn to these. 



V. INTRINSIC DECOHERENCE: THEORETICAL 
FRAMEWORK 

As noted earlier, by 'intrinsic decoherence' we mean a 
process which looks like decoherence (ie., one in which 





FIG. 5: (Color on-line) Contributions to Q{R,R') from differ- 
ent paths. In (a) we show the contributions in conventional 
quantum mechanics - paths sum independently. In (b) we 
show the intrinsic decoherence contribution, in which the cor- 
relator k(1,2) connects different branches of the propagator 



phase coherence or phase interference is destroyed with 
time) but where this process is intrinsic to Nature, ie., 
not obtained by simply averaging over 'external' degrees 
of freedom that happen to be entangled with the system 
of interest. As such, an investigation into intrinsic deco- 
herence amounts to a search for a breakdown of quantum 
mechanics, of a kind that mimics (at least up to a point) 
conventional environmental decoherence. Thus we are 
looking for a breakdown of quantum mechanics, of a spe- 
cific kind. 

For many physicists the first reaction to this idea is 
probably to ask why anyone would bother looking at all! 
The two main reasons are (i) the well-known difficulties 
attending macroscopic superpositions of states in quan- 
tum mechanics [3l| , and (ii) the clash between Quantum 
Mechanics and General Relativity [s^], which has yet to 
be resolved - the success and generality of each of these 
theories makes a resolution of this clash perhaps the most 
important single problem in physics today. 

As noted above, most previous discussions of the possi- 
ble breakdown of quantum mechanics have started from 
the Schrodinger equation, and added extra non-linear 
terms. In what follows I will very bricfiy sketch a differ- 
ent framework. The purpose of this sketch is not to give 
anything like a full exposition of this theory, for which 
there is no space here, but simply to indicate how this 
theoretical framework can be set up. Full details will 
appear elsewhere [ssj . 

Our choice of theoretical framework is motivated by (i) 
the observation that the quantum phase, via its connec- 
tion to the action, plays a fundamental role in the time 
evolution of all known physical systems; and (ii) that the 
transition amplitude (or just the 'quantum amplitude') 
plays a similarly fundamental role. It is hard to sec how 
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these would no longer be incorporated in any new theory 
intended to go beyond quantum mechanics. The empha- 
sis on the quantum phase also focusses attention on phase 
coherence. 

Our starting point will be the standard Feynman path 
integral formulation of quantum mechanics, and in this 
short paper we will discuss only the non-relativistic ver- 
sion of this, for ordinary particles. The great advantage 
of the path integral formulation is that it makes explicit 
the idea that at any particular point on its worldline, a 
particle can choose between all possible paths. We will 
now generalize this idea as follows. We define a propa- 
gator G{R, R') for a quantum system between two space- 

I 



in which the interpath correlation functional k[xi,X2] (a 
fimctional of the two paths xi(t), X2(<)) remains to be de- 
termined. The above definition is heuristic in the same 
way that the usual path integral is, in that the normal- 
ization has yet to be determined (ie., the measure for 
the paths has yet to be defined): this point is discussed 
below. The system is clearly non-relativistic - we have 
separated space and time, and the correlator k[1,2] de- 
pends only on integrals over the spatial variables xi,X2, 
so that the correlations it mediates are instantaneous. 

The correction term AQ{R,R') describes pairs of 
paths, and it allows communication or correlations be- 
tween different paths. The physical interpretation of ([23]) 
is not that that of a self-interaction - in conventional 
quantum mechanics, self-interactions for a particle are 
described by single particle paths, in which an interac- 
tion leaves and then rejoins the particle path - one then 
sums over the different possible paths. Here we arc sum- 
ming over all possible pairs of paths between the initial 
and final state, with the correlator k[1, 2] weighting pairs 
of paths (see Fig. [S|); the total transition amplitude is 
now produced by summing not only over paths but over 
correlations between them. We can if we wish think of 
this as some kind of communication between different 
branches of the wave-function. Clearly there is no par- 
ticular reason for stopping at pairs of paths - indeed it 
would be more elegant to assume a sum over all possible 
n-tuplcs of paths, with < n < iV and N ^ oo. How- 
ever we will not discuss this more general formulation 
any further here, since we are simply going to assume for 
the moment that X[2, 1] is very small, so that we can to 
first approximation ignore all higher 7i-tuples. We then 
write 

k[1,2]=7J[1,2]; (7 « 1) (24) 



time points R' = (r',t') and R = (r, i), which we will 
assume to have the form 

g{R,R') = Go{R,R') + Ag{R,R') (21) 

where Go{v,y' \t^t') is given by the usual path integral 
expression in non-relativistic quantum theory, viz., 

Go(r,r';i,t') = ^ 2?x(t) exp ^ dTL(x,x;T) (22) 

whereas the new term Af/(r, r'; t, t') will be taken, in the 
first instance, to have the form: 



(23) 



where 7 is dimensionless and we note that both J [1,2] 
and k[1,2] have dimensions of inverse length (or inverse 
interval in a relativistic version). We will also typically 
assume that J[l, 2] takes the form 

J[xi,X2] =cxp ix[xi,X2] (25) 

where the "phase" x[l: 2] will in general be complex. 

In this theoretical framework, we will assume that both 
the amplitude V'(r, t) and the density matrix still appear; 
however the normalization of the path integral is not so 
obvious, and in fact has to be handled rather carefully. 
The amplitude obeys the equation of motion 

^(r,t) = j rfr'g(r,r';t,0^(r',0 (26) 

However the density matrix is no longer given by 
p(ri,r2,t) — '0(ri, i)'0^(r2, 0- This is because we must 
allow correlations between different paths of the den- 
sity matrix, in line with the assumptions made above 
for the 1-particle propagator. This means that we can 
still write the equation of motion for the density matrix 
p{X,Y) = p(x,y,t) in the form p{2) = J dl/C(2, l)p(l) 
as before, but now the propagator IC{X,Y; X'Y') con- 
tains correlations between all possible paths connecting 
the end points X,Y, X' ,¥' . If we then make the same 
expansion in powers of A that is implicit in (|23p above, 
we get an expression for the propagator of the form 

IC{X, ¥■ X'Y') = K{X, Y; X'Y') + AK.{X, Y; X'Y') 

(27) 

where K{X,Y;X'Y') = g{X, X')g* {Y,Y') is the prod- 
uct of the two independent propagators (each contain- 
ing internal corrections from the correlator k[1,2], as 
above), and the correction AJC{X,Y; X'Y') contains 



Ag{r,r';t,t') ^ X'xi(t) J X'x2(t) k[xi, X2] exp ^ ^ dr [L(xi,xi;t) + L(x2,X2;t)] 
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cross-correlations between the forward and reverse paths 
in the density matrix propagator, with the lowest-order 



J 



term having the expression 



A/C(X,y;X'r') = / Px(t) / Vy{T) K[x,y]exp- / dr [L(x,x;t) - i(y,y;T)] 

•J ^C' ^ Y^' ^ t' 



(28) 



and so on. This last term, because it correlates the 2 
parts of the density matrix propagator, can then have 
the effect of generating what looks like decoherence. 

One can also define a reduced density matrix in the 
case where one deals with a system for which, for one 
reason or another, it is necessary to average over some of 
the degrees of freedom. Suppose the total density matrix 
(for the 'universe') is p(Q,Q',i), where the coordinates 
Q — (R, {xfe}), with k = 1,2,...A^, and we are only in- 
terested in the 'central coordinate R. Then we simply 
write that the reduced density matrix is given as 

p(Ri,R2,t) = rrfBni,] p(Q, Q',i) 

= n / dxkp{(^M',t) (29) 
k •' 

in the usual way. 

We may then, by analyzing the dynamics of the re- 
duced density matrix, examine the corrections to its dy- 
namics coming from the extra term A/C(2, 1) in the prop- 
agator of the full density matrix. Suppose, for exam- 
ple, that we have a Lagrangian of form >C(R, {z/j}) = 
Lo(R)-l-L/(R, [xk]), where the 'fast' Lagrangian is writ- 
ten 



L/(R, {xk}) = L^{{xk}) + Li„t(R, {.Tfc}) 



(30) 



with the fast variables now playing the role of the en- 
vironment. Then we can also make such a separation in 
the propagator C7(2, 1). The 'bare' term has the standard 



J 



form [34 1 



/•2 

Go(2,l) = y PRe*^''*^°(^'*)G{(2,l) (31) 



where G^(2,l) = Gi{{xX\xX^yMM\Mt)]) is the 
propagator for the fast variables, conditional on some 
path R(t) for the slow 'central' variable (ie., the propaga- 
tor for the {xk{t)} when the slow variable path is 'frozen' 
to be some specific R(t)). Typically one sandwiches the 
fast propagator between the instantaneous eigenstates 
7i(R(t))) of the entire collection of fast variables, and 
then we can write that 

G{(2,1) = (2|n(R))G{(TO,n;R(0)(m(R)|l) (32) 
where G{(n,m;R(t)) = e^k ! ^t' Li^^{-R.{t') ^ ^j^-j^ 
five Lagrangian 

iJ;„(R(i')) = io(R) - e„(R) - iKk ■ (n|VRm) (33) 

containing the usual 'Berry phase' term derived by Born 
and Oppenheimer. 

Now suppose we add the extra term /S.Q{2, 1) to the 
propagator. If we again make the separation between 
slow and fast variables, we get an extra phase term in the 
dynamics of the Green function, and similar extra terms 
in the propagator for the density matrix, whose effect is 
to correlate the variables between 2 separate paths, for 
both the fast and slow variables. Substituting ([50)1 into 
(l23l) . we find that these phase factors have the form 



$„™[R,R'] ^ y"diR- (n(R)|VR|a(R)) Xa;3[R,R'] (a(R')|VRHm(R')) • R' 



(34) 



r 



where Xap [Rj R'] is given by 

Xo,3[R,R'] = |a(R)) '^[R.R'] (/3(R')I 



(35) 



In analyzing a macroscopic system with many degrees of 
freedom, this kind of technique is essential. 

One may also look at various limiting cases in this for- 
malism. The semiclassical limit ?i — can be studied 



directly in (|23|) . and we see that since fi appears in the 
denominator of phase factors, the same subtleties arise as 
for ordinary quantum mechanics. In the same way one 
may study the adiabatic limit, both for a single parti- 
cle and for an A^-particle system, and the 'macroscopic' 
limit, where the number of particles in some object is 
taken to be very large. 
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Up until now I have said nothing about how we are 
supposed to interpret the density matrix, and expres- 
sions Hke {A) = Tr {pA\ (in standard quantum mechan- 
ics of course these refer to expectation values and lead, 
via the Born rule, to probabilities for these). Quite apart 
from these interpretational questions, we see that there 
must be a normalization requirement on the functional 
integrals - this is of course a problem even in the path 
integral formulation of ordinary quantum mechanics. In 
the present case the normalization (and the interpreta- 
tion of the interpath correlation) depend on what specific 
form we employ for the kernel k;(2, 1); this has then to be 
determined for each case. This should not be too surpris- 
ing - indeed, the measure of paths in the path integral 
in ordinary quantum mechanics also depends on which 
specific Lagrangian we choose. To give a full discussion 
of this would take us into too long a technical diversion 
33|. 

At first the modifications of ordinary quantum me- 
chanics seem rather small, if the kernel k(1, 2) is assumed 
to be small. However the effect on the dynamics of any 
quantum system actually depends critically on the form 
of the kernel k(1, 2). It is then interesting to explore one 
example, which is of considerable interest in its own right. 
Thus we now turn to the specific mechanism of gravita- 
tional decoherence, in which a correlation between paths 
mediated by gravitational interactions is assumed to arise 



VI. AN EXAMPLE: GRAVITATIONAL 
DECOHERENCE 

Our example of an intrinsic decoherence mechanism is 
taken from studies in quantum gravity. Models such as 
the one to be discussed have to be taken seriously, given 
the apparent incompatibility of between Quantum Me- 
chanics and General Relativity. However I immediately 
emphasize that what follows is not intended to do more 
than indicate, by way of a toy example, the sort of thing 
that might be done - it is very far from being any kind 
of polished theory. 

The basic physical idea to be discussed is that gravita- 
tional effects must lead inevitably to some sort of intrinsic 
decoherence in the time evolution of isolated systems. Wc 
see that any such arguments lead us inevitably to a non- 
unitary time evolution for state vectors - no averaging 
over an environment is involved, and yet phase correla- 
tions are lost. Such proposals have been quite common 
in the literature, ever since Hawking proposed that such 
non-unitary time evolution ought to be a part of any 
quantum theory of gravity [s^. Some of these propos- 
als have incorporated quantum spacctime fluctuations at 



the Planck energy as the source of this intrinsic deco- 
herence [3^, even to the extent of estimating how these 
might cause apparent decoherence in the dynamics of a 
superconducting SQUID . These theoretical proposals are 



controversial jSTl |38| . and there has certainly been no 



clear experimental test of them at all (note that the esti- 
mates given in the theory for intrinsic decoherence rates 
contain a number of factors ~ 0(1), which then end up 
being exponentiated). 

Somewhat distinct from these 'Planck scale' propos- 
als are the ideas discussed by Penrose [l^], following on 
from Diosi [2lj. These authors argue that in analyzing 
quantum superpositions of different spacetimes, coming 
from situations where a mass is in two different posi- 
tions, a kind of 'gravitational time-energy uncertainty 
principle' must operate. In these analyses, the timescale 
Aic = h/ /S.Eg plays the role of a decoherence time; 
here A£'g is the energy scale associated with the gravi- 
tational self-interaction of an object between its different 
quantum states. However this decoherence is not envi- 
ronmental - it can be viewed as either coming from an 
uncertainty in the background time, or as arising from 
a mismatch of the spacetimes generated by the gravi- 
tational field of the system in question, in its different 
quantum states. 

How can we incorporate these ideas into the formalism 
above? Let us, for the sake of having a concrete example, 
consider a correlator k(1, 2) of form 



K[r,r'] = exp j dr 



|r(T)-r'(r)| 



- 1 (36) 



so that the strength of the phase communication between 
different paths is now proportional to the strength of the 
'gravitational self-interaction' between them. I make no 
attempt here to justify this choice for K[r,r'] except to 
note that it is consistent at least with the spirit of the 
Diosi-Penrose uncertainty arguments above. We have as- 
sumed that the gravitational self-interaction can be han- 
dled in a Newtonian approximation. 

Note, however, one key difference with the simple form 
introduced above - the correlator K[r, r'] is no longer mul- 
tiplied by a small parameter A. Indeed the phase factor 
actually blows up at short distances, even when the mass 
TO is small; and no power series expansion of the expo- 
nential is valid in the short distance limit. 

Consider now a really simple problem, in which we 
have a free particle of mass to propagating from X' = 
(R',i') to X = (R, t). In the absence of the correc- 
tion AQ, wc just have the usual free propagator, ie., 
g{X,X') GoiX,X'). Adding the gravitational term 
then gives the correction 
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Pxi(r) / 2?X2(r) K[xi,X2]exp— / dT—{±l+£ 



(37) 



r 



This correction affects only the relative coordinate r 
Xi — X2 , and thus we see that 



Ag(X, X') oc A{Q, 0; i, t')Go{X, X') 



(38) 



where ^(r,t;r',t') is just the return amplitude (starting 
and finishing at the origin) for a particle of inertial mass 
m in an attractive 'Coulomb' field, with a 'charge' of 



strength ^TiGm? . Thus the main effect on the 1-particle 
propagator is just a finite renormalization, whose effects 
are not obviously important, since we expect them to 
be absorbed in the normalization of the path integral. 
However the effect on the density matrix propagator is 
more dramatic - we have a correction term which we write 
out in full, ie.. 



A/C(X,F;X'F') 



2?x(t) / 2?y(r) 



exp / dr 



|x(r) -y(r) 



- 1 



xexp^ /dT^(x2-y2) (39) 



in which decoherence is now caused by what looks like a 
gravitational interaction between forward and backward 
paths of the density matrix. Indeed, the kernel K[x,y] 
is now acting as an infiuence functional. Decoherence in 
this calculation will appear in the reduced density matrix 
in both the momentum and position representations. 

An obvious question one can ask about this calcula- 
tion is that, as noted above, the intrinsic decoherence 
correction to the dynamics here is not small in the limit 
|x(r) " y(''')| — > 0. Thus we have no reason to believe 
that these calculations, limited as they are to the lowest 
order in the expansion in k[1, 2], give a proper picture of 
what is happening. This might be alright if the corre- 
lator as given above is assumed to be exact, but this is 
hardly likely. Another question, that has to be dealt with 
if one wishes to discuss possible experiments, is to gener- 
alize this calculation to deal with a system of many par- 
ticles. Thus one can imagine doing an experiment with a 
large mass made from particles distributed according to 
some density distribution. The way to deal with this is 
precisely using the separation of slow and fast variables 
discussed above - the calculations are quite lengthy [33| . 

This if course leads us to the key question - how easy 
will it be to see departures from quantum mechanics 
caused by intrinsic decoherence? A key problem with 
many of the theoretical attempts to predict a break- 
down of quantum mechanics is that, because they are 
attempting to describe some kind of wave-function col- 
lapse, they give predictions for dynamics which look very 
much like standard decoherence. Thus one can ask - how 
can we distinguish between intrinsic and environmental 
decoherence in experiments on, eg., a macroscopic system 
like a superconducting SQUID, or an array of magnetic 



molecules? 

This is clearly a research problem of great interest. On 
the experimental side, some very ambitious and inter- 
esting efforts exist to look for gravitational decoherence 



40|, in which the aim is to set up quantum states in 
which reasonably massive objects are in superpositions 
of spatially separated states. One very nice feature of 
these experiments is that the gravitational decoherence 
rate should indeed vary with external parameters in a 
way rather different from what one would reasonably ex- 
pect from environmental decoherence rates (at least those 
coming from currently understood environmental deco- 
herence mechanisms). However the problem, at least for 
tests of gravitational decoherence, is that there seems 
to be no unambiguous way of calculating the gravita- 
tional decoherence rate using the simple time-energy un- 
certainty arguments that the Diosi-Penrose arguments 
use. The answers depend in a quite arbitrary and very 
sensitive way on how one 'coarse-grains' the density dis- 
tribution of the objects involved - this is simply because 
of the divergence of the gravitational 'self-energy' term 
as two mass distributions approach each other, so that 
the behaviour of the density distribution at very small 
distances determines the intrinsic decoherence rate. One 
interest of the above approach is that it may help us to 
clarify this problem. 

It is clearly going to be a very important task for the- 
ory to come up with a clear quantitative prediction for 
gravitational decoherence rates - and in fact since exper- 
iments of this kind are perhaps not too far from comple- 
tion, there may not be too much time left to do this! 
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